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Abstract 



We address the deformation quantization of generally parametrized systems displaying a 
natural time variable. The purpose of this exercise is twofold: first, to illustrate through a ped- 
agogical example the potential of quantum phase space methods in the context of constrained 
systems and particularly of generally covariant systems. Second, to show that a causal repre- 
sentation for quantum phase space quasidistributions can be easily achieved through general 
parametrization. This result is succinctly discussed. 
43 ' 
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^1 Introduction. 

a ■ 

_ _ i 

Generally covariant systems are a particular kind of dynamical systems in which time is included 
among the canonical variables (I] . The state of the system evolves along orbits parametrized by an 
unphysical scalar parameter. The formalism is invariant under reparametrizations of this parame- 
ter (for this reason these systems are also called general parametrized systems) leading to a gauge 
symmetry, which in the Hamiltonian formulation is implemented by a first class (Hamiltonian) con- 
straint. One of the most striking properties of generally covariant theories is that the Hamiltonian 
is identically zero on the constraint hypersurface. Consequently there is no standard Hamiltonian 
time evolution (which coincides exactly with the gauge transformation). Instead, dynamics is to be 
found among the relations between the canonical variables that are determined by the Hamiltonian 
constraint [21 El E] ■ In other words, dynamics is the unfolding of the gauge transformation. The 
most famous example of a generally covariant theory is, of course, general relativity El H] • 

Upon quantization these systems display a number of technical and conceptual problems which, 
motivated by the quest to quantize general relativity, have been intensively studied in the literature. 

* ncdias@mail.telepac.pt 
* joao.prata@ulusofona.pt 
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Since the publication of Dirac's seminal work on the quantization of constrained systems [8], several 
quantization programs have been developed with the aim of refining Dirac's approach and making 
it suitable to address the generally covariant case [U El 13 EH EE] • In spite of this, several impor- 
tant problems concerning the quantization of these systems (the problem of time, the observables' 
problem and the measurement problem, just to mention a few) are still lacking a definitive answer 

EH SEED]. 

One quantization program that has only been scarcely explored in this context is the deformation 
approach J2]-[22]- And this regardless of the fact that the deformation quantization method displays 
some remarkable features that make it especially suitable to address the generally covariant case 
[QUEUES]- On the one hand the deformation approach leads to a formulation of quantum mechanics 
in terms of classical like objects. The theory lives in phase space and mimics the structure of classical 
statistical mechanics. The state of the system is described by a phase space quasi-distribution 
and observables are also functionals on the phase space. From the deformation point of view 
quantization amounts to the substitution of the standard product of functions by a non- commutative 
star product [2UI22]- The theory has been proved useful when addressing a wide range of fields of 
research ranging from topics in non-relativistic quantum mechanics jT^j, [211-121] to some current 
developments in M-theory [3UJ |M1 E2] . In the context of generally covariant systems it has been 
advocated that the relation between classical statistical and standard operator quantum mechanics 
should be emphasized when approaching the quantization of these systems jHS|- The point of view is 
that the problems besetting their quantization, being to a large extend conceptual ones, may receive 
key physical insights from the more intuitive classical statistical analysis. From this perspective the 
deformation methods seem to be especially suited. 

On the technical side the deformation approach acquired the status of a powerful mathematical 
theory not displaying any of the subtleties of the path integral or canonical quantization prescrip- 
tions. Most remarkable is the fact that the theory is able to generate the quantum version of a 
classical system living on a generic (possibly curved) Poisson or sympletic manifold [3U ESI EM EZj ■ 
Of special relevance for generally covariant systems is that the deformation approach displays some 
powerful tools to address the quantum eigenvalue problem, the most meaningful of which being 
a covariant formalism unifying all distinct phase space representations of the eigenvalue problem 
[2UEE1EH], a formal solution of a generic stargenvalue equation [1211101 and a set of efficient methods 
to determine the semiclassical expansion of the stargenf unctions |4"T| I42j . 

In this paper we aim at providing a pedagogical incursion into some of these methods by studying 
the deformation quantization of a special simple kind of generally covariant systems: those obtained 
by parametrizing an originally non-covariant version of the system. The parametrization procedure 
leads to a new formulation of the dynamics living on an enlarged phase space where the physical 
time is incorporated as a canonical variable |1. The existence of such a natural time variable 
considerably simplifies the (technical and conceptual) analysis of the dynamical structure of the 
system. For this reason these models have been used [HEBE] to test a set of quantization techniques 
and interpretative prescriptions aiming at addressing the far more difficult problem of quantizing 
"already parametrized systems" of which the most significant example is general relativity. 

The Hamiltonians of these models (as in the general case) are identically zero on the constraint 
hypersurface. At the quantum level the imposition of this constraint determines the physical states 
of the system. In the deformation context the imposition of first class constraints yields a stargen- 
value equation [H^ |4"3] and so the physical states are the zero stargenf unctions of the Hamiltonian 
symbol. We are then able to apply the powerful tools of the deformation approach to study the 
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quantum eigenvalue problem. A complete characterization of the physical states will be provided 
in three different phase space representations. In particular, we obtain a representation where 
the dynamics of the quasidistribution (with respect to the physical time variable) is dictated by 
the classical Liouville equation. This is an interesting side result that shows that through general 
parametrization a causal formulation for the distributional sector of quantum mechanics is made 
possible. More generally we find that the remarkable parallelism between classical statistical me- 
chanics and phase space quantum mechanics carries on intact to the generally covariant setting, 
a property that establishes a promising starting point for the future research on the deformation 
quantization of generally covariant systems. 

2 Covariant Wigner quantum mechanics 

The aim of this section is to review the main structures of the covariant formulation of the defor- 
mation quantization procedure. Our analysis will be restricted to the case of flat phase space as we 
will be dealing only with systems of this sort. The reader should refer to [361IHE] for a more detailed 
presentation of these results and to [HE1 EZ| for the generalization of the formalism to the non-flat 
case. An important part of this section (eq.(9) to eq.(15)) focuses on the covariant generalization 
of the *-genvalue equation, once again, just for the flat case. This is a topic that will play an 
important role latter on. 

Before proceeding and to avoid future misunderstandings let us make the following remark: 
Throughout the paper the word covariant will be used in two different contexts: to designate 
the invariance under time reparametrizations of the original classical system and to designate the 
invariance under general coordinate transformations of the Moyal plane. The aim of this paper is 
to present the deformation quantization of a generally covariant (parametrized) classical system. 
We use the covariant version of the deformation quantization procedure (a subject where this paper 
makes no original contribution) because this leads to a larger set of possible quantum phase space 
representations, among which we will find the causal representation presented in section 4. 

Let us then settle down the preliminaries: we consider a iV-dimensional dynamical system, its 
classical formulation living on the flat phase space T*M. A set of global canonical coordinates 
{{qi,Pi,i = ^--N}) can then be defined on T*M in terms of which the sympletic structure reads 
w = dqiAdpi. Upon quantization the set {g^} constitutes a complete set of commuting observables. 
Let us introduce the vector notation q = (g 1; q N ) and designate by \q > the general eigenstate of 
q associated to the array of eigenvalues q and spanning the Hilbert space 7i of the system. Let also 
A(H) be the algebra of quantum observables over Tt and A{T*M) the algebra of classical functions 
over the classical phase space T*M. Let (Q, P) be a second set of phase space coordinates, related 
to (q, p) by a generic phase space diffeomorphism (not necessarily canonical): q = q(Q,P) and 
P = p{Qi P)- The covariant Weyl transform [3*Sj : 

A(T*M); (1) 
W^ P) (A) = h N J d N x J d N ye- i ^ p >y5(x-q(Q,P)) < x + ^y\A\x - \y > 

where \x ± \y > are eigenstates of q, yields the entire structure of covariant phase space quantum 
mechanics. The map W^ffl can be applied both to an observable A as well as to the density 
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matrix \ip{t) >< ip(t)\. In the first case it yields the %%-Weyl symbol A'(Q, P) = W^%AA) of the 
original quantum operator and, in the second case, the celebrated Wigner function fw(Q, P,t) = 
(2^F W (|5)(I^(*) >< ^(*)D- Notice that if (Q,P) = (q,p) then the covariant map W^l reduces 
to the standard Weyl transform in the variables {q,p)- In this case we will use the notation Wi^a 
to designate W$$. Notice also that A'(Q,P) = A(q(Q , P) , p(Q , P)) where A(q,p) = W m (A) 
and the same relation is valid for the Wigner function. The covariant Weyl map implements the 
transformation (q,p) — > (Q,P) as a coordinate transformation in quantum phase space and defines 
a covariant star-product and Moyal bracket [, ]m'^ _ EUEBj (A, B e A{7i)): 

which display the functional form: 

A'(Q,P)*' aP) B'(Q,P) = A'(Q,P)e^<ln V] B'(Q,P) 

A'(Q,P),B'(Q,P)] mI _ = y(Q,P) S m(^v\jg P) V^jB'(Q,P), (3) 

where the covariant derivative is given by (let O n = Pi,O l = Pi,i = 1, • • • , N; H = Qi-N, O l = 
qi _ N ,i = N + l,---,2N): 

^A' = & t A', Vy j A'^d f i d' j A'-rp k A', d'^d/30' 1 - i,j,k = l,---,2N, (4) 

and 

j'v (o p) - ^i^ijM j* (d p) - d ° n d2 ° b (5) 

j (q,p)^' ' dO k dO l m ' L ^' r) dO h dO'WO' k ' { > 

are the new symplectic structure and Poisson connection associated with the coordinates (Q,P), 
respectively Finally, notice that in eq.(5) we explicitly took into account the phase space flat 
structure. 

When formulated in terms of these structures Wigner mechanics becomes fully invariant under 
the action of general phase space diffeomorphisms: 

A'{Q, P) *' m B'(Q,P) = A(q[Q,P),p(Q,P)) * m B(q(Q, P),p(Q,P)) V A ,BeA(T*M), (6) 

the covariant generalization of the Moyal and stargenvalue equations reading: 

p' a ., b (Q,P)*'^ P) A'(Q,P) = bp' a;b (Q,P) 

where p' a . b (Q,P) is the a-left and 6-right -^-genf unction of A'(Q,P). These equations transform 
covariantly under arbitrary phase space diffeomorphisms yielding, in any coordinates, identical 
mathematical solutions and thus identical physical predictions: 

P(A'(Q,P;t) =a)= J d N Q J d N P(detjg p ^p' a . a (Q,P)f{ v (Q,P-,t). (8) 



Let us consider the stargenvalue equation (7) in more detail. Let {Ax, A^} and {Bi, -Bjv} 
be two sets of commuting observables satisfying = ihSij. Let A' { (Q, P) = W^flAAi) and 

B[{Q,P) = W^fpABi). Then the simultaneous solution of: 



pL.,a N M,..M(Q^)*{Q,P } ^(Q^P) = b iPai ,.., aN ., h ,.., bN (Q,P) , i = l..N 

is given by [TH] : 



(9) 
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ffai t ..,a N ;b li ..,b N (Q,P) = ■■;A' N ;a 1 , ...,a N ;b x , ...,b N ) = Yl* r p'^.^Q , P) , (10) 

{Q,P) . (Q,P) 

1 = 1 

„i vM {Q,P) s - iJ^^ I (ID 

(Q,P) ATT (Q,P) w ' ' J cq.p) 



where TT^i *' is the iV-fold starproduct and the starexponential e ** ii *^ ,p ^ = Ej (k, Q, P) is de- 
fined as the solution of the differential problem: 4tE*' (k, Q, P) = %A\ (Q, P) ». E*> (k, Q, P) 

" K (Q,P) \Q°! (Q,P) 

and E*i _ ^ (0, Q, P) = 1, which is unique for the Weyl symbol A[(Q, P) of a generic self-adjoint op- 

(Q,P) 

erator A4 (a result that is just the Weyl-Wigner translation of an equivalent result in operator 
quantum mechanics jSj). Furthermore _ ^ (k,Q,P) = J2t=i ^rH^Li*' - - A*i{Q)P) whenever k 

(Q,P) (Q,P) 

is in the radius of convergence of the series, which justifies the notation used in eq.(ll). We also 
introduced the *-delta notation A*/ (A\, A' N \ ai, .... a N \ 6 1? b N ) that will be extensively used 

(Q,P) V 

later on. Notice that the stargenfunctions Pa 1) .. )a ^ ; 6 1) .. ) 6 w (Q, P) are the |^^.-Weyl transform of the 
projectors \a±, ..,a N >< 6 1; ..,6jv| where \a>i, ..,a N > and ..,b N > are simultaneous eigenvectors 
of Ai, i = 1..N. Moreover the stargenfunctions p' a .. b .(Q, P) are the ^jfL -Weyl transform of the 
projector: 

>< 6i| = / dz\a u z>< h,z\ = — [ dke^~ a ^' e ik{A ^\ (12) 
J 2n J 

where z — (oi, Oi-i, Oi+i, ■■■,cin) is the iV — 1 array of degeneracy indices. The Weyl transform 
of the first identity in eq.(12) yields the inverse relation of eq.(10): 

Pa^Q'P) = j rfa i--^-i rfa m--^«ivpl 1 ,.., a! ,.., ajv;ai ,.. ifeli .., ajv (Q, J P)- (13) 
Finally, the diagonal elements p' ai - ai (Q,P) are given by: 

p' ai ; ai (Q,P) = ^ f dkef^ P) ' ai) = A, (4; 0*00 = A, (4- a,), (14) 

LlX J (Q,p) (Q,p) (Q,P) 

this being the object that enters the probability functional (8). The stargenfunctions (10,11) trans- 

— * — * 

form as scalars under arbitrary phase space coordinate transformations. For instance p' a .. a .(Q,P) 
satisfies: 

p' a ,^(Q,P) = A» ; (4(Q, P) - = A {Ai{q{Q,P),p{Q,P)) ~ a,) = Paim {<I{Q , ?), jt($ , ?)) 

(15) 
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where p a .. a .(q,p) is the diagonal solution of the stargenvalue equation in the (q, ^-representation. 

Lastly, we should point out that the former results (in the presented form) are valid whenever 
we are able to provide the complete sets of commuting observables {A{\ and {Bi}. In fact this 
restriction can be considerably weakened. It is crucial to the overall approach that the complete set 
{Ai} exists, but the requirement on the existence of the set {Bi} can be discarded, while preserving 
the validity of eqs.(8,10,13,14) exactly and that of eqs.(ll,12) under slight modifications. This is 
the case, for instance, when Aj displays a discrete spectrum. The reader should refer to [19J for a 
detailed presentation of these results. 



3 Deformation quantization of a parametrized non-relativistic 

system 

Let us consider an arbitrary N- dimensional dynamical system with configuration variables q = 
(qi, (/at), described by the Lagrangian L (q, -£). Starting from the standard action S = J dtL (q, 
we impose the time reparametrization invariance by introducing an unphysical time r and promoting 
t to a configuration variable. The action is now re- written as PQ: 

S = fdTiLo(q,fj (16) 

where the dot represents the derivative with respect to r. The Legendre transform yields the 
Hamiltonian formulation living in a (2N + 2)-dimensional phase space spanned by the canonical 

variables t } P t ,q,p where P t = dL/di, pi = dL/dq\, i = 1..N and L(t } i } q,£j) = tL (q, j). This, as 
expected, is a "zero-Hamiltonian" system: 

H = X<P , <f) = P t + H (q,p) (17) 

where is a first class constraint, A is a Lagrange multiplier and H is the Hamiltonian associated 
to the original Lagrangian L . 

Dirac's quantization procedure imposes the constraint as a restriction on the space of physical 
states: <p\ip >= 0. In the density matrix formulation this equation reads: 

01^x^1 = 1^x^10 = 0. (18) 

To find the complete specification of its solutions, and according to the general method of section 
2, the first step is to introduce the complete set of observables {0, A = [Ai, ...,A N )} and the set of 

operators B = (B 1 , B N ), generators of translations in the spectrum of A, satisfying [Aj,B k ] = 
ih5jk and [0, Bk] = 0, j, k = 1..N. We easily find that: 

Aj = Fjfaft-t) , B k = G k (qJ,-t) (19) 

where Fj(q,p,t) and Gk(q,p,t) are the solutions of the original Heisenberg equations: 

r)F BC 1 

-^ = [F V H ] , -^ = [G k} H ] , j,k = l..N (20) 
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together with the initial conditions Fj(q,p,0) = cjj and Gt(q,p,0) = pk, satisfy the aforementioned 

— * — * 

requirements. A and B are the quantum histories of the system. 

The phase space representation of eq.(18) is determined by the Weyl transform. We may write: 

*(t,<t>,A,B) fw(t, 0, A, B) = f w (t, 0, A, B) = (21) 

by using the map W^ t ^ Xb) or equivalently: 

0(t, P t , q, p) *( t ,p t ,q,p) fw(t, P t , q,p) = fw(t, Pt, q, P) *(t,p t ,q,f) Pt, q,p) = (22) 
by using the map W( t ,p t ,q,p)- We start by considering the simplest representation, which is provided 

— # — # 

by eq.(21). In this case the Wigner function is a left and right zero (t, 0, A, £?)-stargenfunction of the 
constraint symbol 0. To determine fw explicitly we will follow the procedure described in section 2. 
From eqs.(10,ll,14), the fundamental zero stargenf unctions of the Hamiltonian constraint (which 
will be designated by p d ^j are given by: 

p, rb (t, 0, A, B) = A w>a) (0, A; h = 0, a; h = 0, b) = A w> ^ (0) A W>J) (A; a; b) = 

N [ ' b- B 1 

j = l 



N 

S W*(t*A3) II {e*<*->)** * ^ 5(A, - b 3 )} 

3=1 

N ( y +°° I 



-ih 



,=il27r^n! V 2 

5(0) - 



/ dke ik{A i- b * 



N 



3=1 
N 

w n 



i 



5(0) J] <| e^- a ^^- y dJfee**^- 6 '"^-^)} 



a + b, 



(23) 



where lT^=i* (t 0l 5) stands for the iV-fold starproduct and Ilf =1 for the standard iV-fold product of 
functions. The stargenf unctions p Sb -(t, 0, A, B) are simultaneously solutions of eq.(21) and: 



A? *(t,<t>,A,B) 



p srb (t,<f),A,B) = a jPS)$ (t,(f>,A,B) 
P S jltt,<l>,A,B)* MS) Aj = b jPSi$ (t,<l>,A,B) ,j = l..N 



(24) 



Notice that the diagonal elements a = b fully identify a history of the system: pd,a(t, 0, A, B) = 
5(0) UjLi S(Aj — ctj). It is also straightforward to realize that the fundamental stargenf unctions (23) 
are the (t, 0, A, B)-Wey\ transform of the projectors: 



N 

\h = 0,axh = 0,b\= A(0) J] {eWi-^AiAj-bj)} 



(25) 



3=1 



where the general ket satisfies: 

<f>\h, x >= h\h, x > and Aj\h, x >= Xj\h, x > (26) 

and 

A(0) = J dx \h = 0, x >< h = 0, x | = y J dke ik *, (27) 

is the operator analogue of (14). The most general solution of eq.(21) is a linear combination of the 
fundamental solutions (23): 

f w (t, 0, A, B) = — J da J dbC(a)C*(b)p s fi, 0, A, B) (28) 

where C(a) obeys to the normalization condition that is induced by the normalization of the Wigner 
function. Let us then calculate its norm: 

J J J J dtd(f)dAdBf w (t,(f),A,B) 

= jdnw III dtdAdS S dS S Sc(3)c *^ ri' r: '' '" /; ^ - ^y^) = 

= JdtJdaJ dbC(a)C*(b)5(b- a) = J dtj daC{a)C*{a) (29) 

The divergence of the previous integral indicates that we are integrating over the gauge orbits thus 
spoiling the normalization of the Wigner function. In this case this is quite simple to correct. In fact 

— * — * 

we just have to introduce the phase space measure d\i = dcpdAdB and use it from now on whenever 
we have to integrate the Wigner function. The procedure corresponds to cutting the gauge orbits 
through a single time hypersurface. With this measure the Wigner function satisfies: 



Jdfif w = l ,Vt (30) 



provided the parameters C(a) satisfy J da\C(a)\ 2 = 1. The proper normalization of the Wigner 
function determines a new phase space measure and a restriction on the factors C(a). These 
parameters display a natural physical interpretation. To see this explicitly let us calculate the 
probabilities for the output of a measurement of A. The general x — (x± , .. ., Xn )-left and -right 
stargenfunction of A (and simultaneously y-left and -right stargenfunction of the constraint) is: 

N 

p y ,x;y,x = A, (0, A; y, x ; y, x ) = A* (<j>-y)H. A. (Aj-Xj) = 

N 

= 5(<f>-y)l[5(A j -x j ) = 5(<!>-y)5(A-x) (31) 

3=1 

Therefore, the probability density for A = x and = y is given by: 

V{<f> = y,A = x) = j dnf w {t,<j>,A,B)5{<f>-y)5{A-Z) = J dBf w (t,y,x,B) 



(2nh) 



N 



dB J J dadbC(a)C*(b)5(y) f[ j eW*- a ' )B '8{xj - 



— * 

a + b s r , , , . <N 2 



dadbC(a)C*(b)5(b-a)5(x — )5(y) = \C(x)\ 2 5(y) (32) 



from where it follows: 

P((f> = y,A = x) = lim f V+e dy'V{<t> = y',A = x) = \C(x) \%, . (33) 

It is clear from the previous equation that the term |C(:c)| 2 represents the probability for the system 

— # — # 

to be found in the history A = x. Under the measurement of A with output x the Wigner function 
will indeed collapse to the state: 

Po,S;o,s^S(<f))5(A-x) (34) 

Let us also point out that in the (t, <fr, A, B) -represent at ion the Wigner function is static both with 
respect to the external time r and to the physical time t. In fact from eqs. (23,28) we see that 
§fw = 0. 

The former results can now be easily translated to the (t,P t ,q,p) representation where the 
intention is to solve eq.(22). The two representations are related by the unitary transformation: 

*(*, P t, &P) = U' 1 *(t,P t ,q,p) t *(t,P t ,q,p) U = t 

4>{t, P t , q,p) = W (t ,P t ,q,p)((f>) = U' 1 *(t,P t ,q,p) Pt *(t,P t ,q,p) U = P t + H (q, p) 

Aj(t,q,p) = W ( t t p tm (A j ) = U- 1 * (jtt i %m q j * ( t t p tm U , j = l..N 
Bj{t,q,p) = W(t,p t ,^(B j ) = U~ 1 *(t,p t ,^Pj*(t,P t m U ' j = l -- N ( 35 ) 
where U = e* ( f^^*.The fundamental stargenfunctions are: 

p s j(t, P t , q,p) = \ ttPt ,J(f), A; h = 0, a; h = 0, b) (36) 

N r i 

J'=l L } 

which in general do not simplify any further as in equation (23). These stargenfunctions are the 
(t, P t ,q,p)-Wey\ transform of the projector (25). The most general solution of (22) is then: 

f w (t, P t , q,p) = -^-^ J dadbC(a)C*(b)p^(t, P t , q,p) (37) 

and can be obtained directly from the Wigner function (28) by applying the unitary transformation 
(35). Hence, fw(t, P t ,q,p) is properly normalized for the phase space measure d/j, = dP t dqdp 
provided the coefficients C{a) satisfy the normalization / daC(a)C*(a) = 1. Notice that in general 
the new Wigner function does not have support exclusively on the classical constraint hypersurface. 
This is due to the non-local character of the *-delta functions in eq.(36). Notice also that fw does 
not evolve with respect to the external time r but it displays the standard time evolution with 
respect to the canonical time t: 

df w (t,P t ,q,p) = [m p ^ ^ f ^ p ^ mM{t p ^ = o 

<=► d J^lMlPl = [H (q,p)J w (t, P t , m\M lM (38) 

We conclude that the (t, P t , q, ^-representation yields the extended phase space Schrodinger picture 
for the quantum generally covariant system. Likewise the (t, 0, A, B) -represent at ion provides the 
phase space Heisenberg picture. 
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4 Causal representation 

We now study another possible phase space representation of the system. The first step is to 
specify a new set of classical phase space coordinates (t, P t ,Q, P). Let us define the phase space 
diffeomorphism by: 

t = t , P t = 4>-H {A,B) , Q = Q(t,A,B) , P = P(t,A,B) (39) 
where H (A, B) = W {Afi) {H Q ) = W m (H )\. =A ^ =s and Q(t,A,B), P(t,A,B) satisfy: 

^ B) = {Q(t, A, B), H (A, B)} ( a,b) , dPit d f B) = {P(t,A,B),H (A,B)} {AS) (40) 
together with the initial conditions: Q(0, A, B) = A and P(0, A, B) = B. That is Q(t,A,B) and 

— * — * — * 

P(t, A, B) constitute the classical time evolution of the deparametrized system. The transformation 
(39) is canonical and can be easily inverted. It yields: 

t = t , <f> = <f>'(t,P t ,Q,P)=P t + H (Q,P) , A = A'(t,Q,P) , B = B'(t,Q,P) (41) 

where this time H (Q,P) = W m (H )\^ Af=p and A'(t,Q,P), B'(t,Q,P) satisfy: 

dAI(t d f P) = {H (Q, P), A'(t, Q, P)} aP) , mt £' P) = {H (Q, P), B'(t, Q, P)} (m (42) 

together with the initial conditions: A'(0,Q, P) = Q and B'(0,Q,P) = P. This means that the 
functions A'(t,Q,P) and B'(t,Q,P) are the classical histories of the system. We should point 
out that in general the two functions A'(t,Q,P) and A(t,q,p) given by eqs. (41,35), respectively, 
as well as B'(t,Q,P) and B(t,q,p), display different functional forms (they are respectively the 
classical and the quantum Weyl-Wigner histories of the system). To make this point explicit, we 

introduced the prime notation which is in agreement with the fact that A'(t, Q, P) = wf '^'^l^A), 

B'(t,Q,P) = W^ S ^(i) and <f>'(t, P t ,Q, P) = W^ S ^(4>). Finally, let us concisely designate 

the transformation (41) by f : R( 2iV+2 ) — ► R( 2JV+2 ); (t, P u Q, P) — ► (t, 0, A, B) = f(t, P t , Q, P). 

The transformation T is canonical but in general it does not preserve neither the starproduct 

nor the Moyal bracket. Using the generalized Weyl transform W^'f,'^ 8 ! we find a new phase space 

{t,Pt,Q,p) 

representation of the eigenvalue equation (18): 

<j>% P t , Q, P) *{ mAP) fw(t, Pt, Q, P) = f{y(t, P t , Q, P) *\ t ,P t ,Q,P) P *> Q,P)=0 (43) 

where the *| t Pj ^ ^ is the covariant starproduct given by (3) with: 

f)n>i f) 2 n b 

j'v _ _ = pJ _ an d r'\ = _ (aa) 

where O b E {t, 0, A, B} and H G {t, P t , Q, P}. The solutions of eq.(43) can be read from eq.(23,28): 



(45) 
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where: 



p' si {t, P t , Q, P) = A #; U'(t, P t , Q, P), A'(t, Q, P); h = 0, a; /i = 0, 6) 

' (t,Pt,Q,P) 



N 

= A*/ _(^,P t ,Q,P))n*; ,A, ; (4(*,^, J P),a J -,6 J -) = 

(t,p t ,Q,p) A , (t,p t ,Q,p) (t,p t ,Q,p) x ■> 

AT 

— * — * 

= A 



*(t,0, A, S) 

j = l 



AT 

n 



= 6(<f>>(t, P t , Q, P))ei$- s > S '^8(A'(t, Q, P) - *±*) (46) 

The new Wigner function satisfies: 

f^t, P t , Q, P) = f w (t, <f>'(t, P t , Q, P), A'(t, Q, P), B'(t, Q, P)) (47) 

where fw is given by eq.(28). Furthermore and since the coordinate transformation is canonical we 
have: 

dp = dcpdAdB = 5(t)dtd(pdAdB = 5{t)dtdP t dQdP = dP t dQdP, (48) 

and the new Wigner function obeys the proper normalization: / dpf w = 1. Also notice that 
both the diagonal and the non-diagonal stargenf unctions (and thus also the Wigner function) have 
support only on the classical constraint hypersurface a property that it is not shared by the previous 
representation (36,37). 

Finally, let us briefly elaborate on the dynamical structure of the system in this representation. 
It is clear from eq.(43) that the Wigner function is static: 

■j^flv(t,P t ,Q,P;T) = [H\t,P t ,Q,P)J{ v (t,P t ,Q,P;r)} M > tp ^ p) 

= A[0'(t, P t , Q, P), f^t, P t , Q, P; t)] m; = 0, (49) 

confirming the typical picture of frozen dynamics. However, the Wigner function does evolve with 
respect to the physical time t. From eqs. (42,47) we have: 

J>,p„Q,p) = ^(f( t , fl ,ft^f + ^(f( t ,i.,, ( J,#))M + g^ ( f (( ,, ! , ( j,^ 

N or N or 

= E iHr(T(t, P t , Q, P)){H , A'} {mAP) + £ J^(f(t, P t , Q, P)){H , B$} mM 

j=l urL 3 j=l UJD 3 

= i H O,fw}(t,P t ,0,P)> ( 50 ) 

reproducing the classical Liouville equation {H' , fw}^ tPt q p) = 0- 



(t,<t>,A,B) 
{t,Pt,Q,P 

picture: the Hamiltonian vector field lives on the extended phase space and is given by 



We see that the )' p I i-representation of the system leads to an interesting mathematical 



d N dH d N dH d 

& = ^« + SWa%"§a^fl^'" (51) 
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The flows of this vector field define lines in phase space. These lines cross each time hypersurface 
once and only once: they are the histories of the system. These histories are one-dimensional 
hypersurfaces that can be identified by the values of the 1 + 2N constants of motion: 0, A, B. 
Along these lines the Wigner distribution function is constant. However, the correlations between 
the canonical variables, namely between Q and t, do change and time evolution is imprint on these 
correlations. The interesting point is that this picture is not of the classical description but of the 
quantum mechanical instead. 

A natural question is then: what happened to the quantum content of the theory? What 
happened to the interfering trajectories and to the non-local behavior? The answer is that the 
quantum features have been completely removed from the distributional sector and are now exclu- 
sively placed in the observables' sector of the theory. This can be checked explicitly by applying 
the proper generalized Weyl transform to one of the fundamental operators of the system. We have 
for instance: 

Kiffp)^ = gi(t ' *> § )\{t,<t>AS)=?imAP) = M> P " & P), A'(t, Q, P), B\t, Q, P)), 

(52) 

leading to the time evolution: 

dt dt d(j) dt pidAj dt ^{dBj dt 

N Q N Q 

= [q u H } M - H o}(Q,P) ~ 12 bW^' H °}(Q> p ) 

j=i u - n -j j=i ulJ j 

= [q^HolM^ - {qi,H } {JtS) , (53) 

where the identity -B- = [qi,H Q ]M, x - follows from eq.(20) (or alternatively from eq.(35)), we used 

the fact that <ft' is time independent and that the transformation (Q, P) — > (A, B) is canonical. 
Equation (53) does in fact constitute a quantum correction to the classical statistical description 
where we have: 

— * — * — * — * 

Notice that in the classical description Q = q, P = p and q(t, (fi'(t, P t , q,p), A'(t, q,p), B'(t, q,p)) = q. 

We finally conclude by pointing out that a causal representation can also be obtained for de- 
parametrized systems by using an explicit "time dependent phase space representation" where the 
starproduct and Moyal bracket are themselves (scalar) time dependent jSH]- in this approach we 
apply a time dependent generalization of the Weyl map to the density matrix formulation of the de- 
parametrized system and obtain a phase space causal representation of the Wigner function, which 
is formally identical to the one described by eqs. (45,46,52), although t is a canonical variable in the 
former equations and an external scalar parameter in the approach of 39j. Accordingly, the two 
resulting quasidistributions live on different phase spaces. In spite of this the two formulations are 
consistent with each other leading to the single conclusion that a phase space formulation of quan- 
tum mechanics where the distributional sector displays a fully classical causal structure is made 
possible through a suitable choice of representation. This may either be a (scalar) time depen- 
dent representation in the deparametrized phase space or a history representation in the generally 
covariant setting. 
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5 Example 



To illustrate our previous results let us consider the simple system compose by two coupled particles 
and described by the Hamiltonian: 

where are the canonical variables of the particle of mass M, (52, P2) those of the particle of 

mass m and A; is a coupling constant. 

The generally covariant version of this system is obtained by promoting t to a canonical variable 
and imposing the time reparametrization invariance. The extended Hamiltonian formulation of the 
system lives on a 6-dimensional phase space spanned by the canonical variables t, P t , q = (q±, q2),P = 
{.V11V2) which satisfy the commutation relations {t, P t } = 1, {qi,Pj} = o~ij, i,j = 1, 2, all others being 
zero. Upon quantization the system is described by the "zero Hamiltonian": 

"2 "2 

H = \^ , ^ = p t + H = P t + ^ + ^ + kq i f 2 (56) 

where is the first class Hamiltonian constraint and A is a Lagrange multiplier. In standard 
the Dirac formulation the physical states of this system are the wave functions ip solutions of the 
constraint equation: 

4>\if; >= (57) 

We now address the deformation quantization of the system. As in the main text, three distinct 
quantum phase space representations will be presented. 



1. The map W( t ,p t ,q,p) an d the generally covariant Schrddinger picture. 
The Weyl map W(t,p u q,p) yields the (t, P t , g, p)-constraint symbol: 

= W {mm (j>) = P t + Jl + || + kq lP j 



(58) 



and the quantum phase space version of the Hamiltonian constraint (57) is given by: 

(f)(t, P t , q, p) *(t,p t ,q,f) fw(t, P t , q,p) = fw(t, Pt, q, P) *(t,p t m 0(*> P *> 9, P) = ( 59 ) 

The solutions of this right and left stargenvalue equation are given by eqs. (36,37) and they do 
not easily simplify any further. It is worth noticing that due to the non-local character of the 
♦-product the phase space quasidistribution fw will have support on phase-space regions which are 
not classically allowed. One can easily calculate the evolution of fw with respect to the physical 
time t. From eq.(59) we have: 



w 



dr 

dfw 
dt 



A[0, fw]M Mim 

{H , fw}(q,p) + 



= 
24 



w 



dt 

2 ^2kp 2 , 



d 2 h 



dq 2 dpi 



[H Q , fw]M( m 

-{2k qi , 



w 



d 2 f w 
dql 



(9,P). 



(60) 
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which is obviously not of the form of the Liouville equation. Consequently, the Wigner function 
does not display a classical causal structure. 



2. The map W^^^g^ and the generally covariant Heisenberg picture. 

Following the approach of section 3 our first step is to determine the quantum histories associated 
to H. We easily find that [check eqs. (19,20)]: 

' M = & -ft* 
B 2 = p 2 

satisfy the requisites following eq.(18) i.e. Aj, Bj, j = 1,2 commute with the constraint and 
with i and furthermore they satisfy the Heisenberg algebra [Ai,.Bi] = [A 2 ,S 2 ] = ih, all other 

commutators being zero. Since [t,(j>] = ih the set {i, 0, A — (Ai, A 2 ), B = (Bi,B 2 )} is a complete 
set of fundamental operators for this system. 

Using the Weyl map W^ t ^ ^ ^ we find the Heisenberg phase space representation of the constraint 
equation (57): 

<t> *(t,t,AB) fw(t, 0, A B) = f w (t, 0, A, B) * MS) = (62) 

— * 

The fundamental solutions of this equation p-j (where a = (01,02) and b = (61,62)), and also 
satisfying: 

A 3*(t,4>,A,B)Pa,b = a jPa,b > Pa,b *(t,4>,A,B) A 3 = b jPa,b ' 3 = h 2 (63) 

are given by: 

Psrb (t, 0, A, B) = S^^^-^^A, - a -^)5{A 2 - *L±k) (64) 

And the Wigner function, solution of (62), is just a hermitian combination of the fundamental 
solutions. We have: 

f w (t, ( j ) ,A,B) = 1 ^-JdaJ dbC(a)C*(b)p^(t, 0, A, B) (65) 

where C (a) obeys to the normalization condition that is induced by the normalization of the Wigner 
function: / da\C(a)\ 2 = 1. The Wigner function satisfies eq.(62) which implies that = 0. We 
also have from eq.(64) that = 0. That is, the Wigner function is static both with respect to 
the unphysical scalar parameter as well as to the physical time. This is an expected result since, 
in this representation, the Wigner quasidistribution is an exclusive function of observables of the 
system (the histories) i.e. of quantities that commute with the Hamiltonian constraint. 

On the other hand, in this representation, the stargenf unctions of the fundamental variables t, q 
and p do evolve with respect to the physical time. For instance (let \x > be the general eigenket of 
qi with associated eigenvalue x): 

g x (t,cf>,A,B) = ^ (t) ^ i5) (| a ;><a;|) = A WjS) (g 1 (t,^,S)- a ;) (66) 

= 77- / dkei kMt _' A ' S) - x) = / dke lkMtX§) - x) = 5( qi (t,A,B) - x) 
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where qi(t,A,B) can be easily derived from eq.(61): qi(t,A,B) = W, t ^ £g\(qi(i, A, B)) = A 1 + 
jfct — j^B^t 2 . Hence, g x (t,(f>,A,B) satisfies: 

d ;9x (t, 0, A, B) = [g x (t, 0, A, B), H (A, B)] M{13) = {g x (t, 0, A, B), H (A, B)} (IS) , (67) 



df 



B? 



where H (A, B) = ^ (H ) = + + kA\B\. We conclude that the Weyl transform ^Xb) 
casts the phase space dynamics in the Heisenberg picture. Accordingly, the time dependence is ex- 
clusively displayed by the observable (stargenfunction) sector of the theory. 



3. The map W^'p^^} and the generally covariant causal picture. 

Following the prescription of section 4 let us introduce a new set of phase space coordinates 
(t, P t , Q = (Qi, Q 2 ), P = (Pi, P 2 )) defined by eqs.(39,40). We get: t = t,P t = <f>- H (A, B) and, 



Qi 
Pi 
Q2 



A i Bi - k p2<2 

Bi - kB\t 

A 2 + {f + 2kA 1 B 2 }t 

, k r> r> + 2 k 2 r 3j.3 

B 2 



A 1 
B 1 

A 2 

Bo 



Pi + kP 2 H 

Q 2 -{% + 2kQ 1 P 2 }t 

, k p p + 2 , A; 2 p3 + 3 
P 2 



(68) 



We notice that Q(t, A, B) and P(t, A, B) coincide with W^ t ^^(5), W^ t ^ Xb)(p)^ respectively (this 
is an easy result that follows from eq.(61)), i.e the classical and the quantum histories of this 
system are exactly the same. Indeed, eq.(68) solves both the Moyal and the Hamiltonian equations 
of motion. In the notation of section 4 we have A' = A and B' = B (one should notice that this 
result is not valid in general). Hence, for this system, we are not required to introduce a second set 
of "classical coordinates" and can simplify the notation by making: Qj = qj and Pj — pj, j — 1, 2. 

We now consider the action of the generalized Weyl map W^'p^^}. The associated covari- 
ant starproduct *'(t,p t ,q,p) and Moyal bracket [, }m { 
pi, O' 3 = p 2 , O' 4 = t, O' 5 = q u O' 6 = q 2 , O 1 = 0, 
and i, j = 1..6) [check eq.(5)]: 
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(t,P t ,q,p) 



fij 



J (q,p) 



are characterized by (let O = Pt, O 
Si, O 3 = B 2 , O 4 = t, O 5 = Ai, O e = A 2 

(69) 



p/i 
1 22 



p/5 

24 



p/5 

42 



J_ p/1 

M ' 1 24 



p/1 
1 42 



p/1 
1 34 



- r' 1 
- 1 43 

p/i 

p/5 _ 
1 34 ~ 



_p/6 
1 44 



p/5 _ 
1 43 ~~ 



2fr 

-MP1P2, 



p/1 

1 35 



p/1 
1 53 



_p/5 _ 
1 44 " 

1 34 



- r 2 ■ 
■ 1 43 ■ 



1 33 _ 

_p/6 
1 45 



2kA 1 (t,q,p), 
2kp 2 , 



p/2 
1 33 



p/6 

35 



p/6 
~ L 53 



all other Christoffel symbols being zero. 

In this representation the constraint equation reads: 



2kt, 

p/6 _ 
1 34 " 



p/5 
1 33 

p/6 . 
1 43 ■ 



p/6 
1 23 



p/6 
1 32 



p/6 _ 
~ L 54 — 

k_ + 2 



-±-2kA 1 (t,$,fi 



(P t + ff (?,p)) *' m m fw& p ti = P *> ^ "W,a ( P * + #o(?,p)) = (70) 

where H (q,p) = W^'p'^^(H ) is given by eq.(54). The fundamental solutions of eq.(70) are given 
by: 

Psfa P *' 9,P) = A *{ t . Pt ^ p t, A (t, q,p);h = 0,a;h = 0,b) 
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= 5Mt, P t , q,p))e^-^^6(A(t, q,p) - *±*). 

= S(P t + Ho(^)e^^^^^5 (,i - - ^ 2 - ^) 

.,( fc -{^ + 2^}* + ^ + ^#-^) (71) 
and the Wigner function is once again a hermitian combination of the fundamental solutions p'_ r : 

fw(t, Pt, Q,P) = J^K^i^ ><* = / rfCWCrtif/^t, P t , g,p), (72) 

It is related with the (t, (ft, A, .^-representation by: 

fUt, P t , q,p) = f w (0, <P(t, P t , q,p), A(t, q,p), B(t, q,p)) (73) 

where fw is given by eqs. (64,65). We conclude that: 1) The support of f' w is confined to the 
classically allowed regions [check eq.(71)] and 2) its evolution with respect to the physical time 
satisfies the Liouville equation: 

d fw _ y- df w dAj 2 dfw dBj _ 2 df w . A dfw r „ R , _ f „ f , x 

~dr-^ 1 dA~^r + ^ 1 dB~^f- £ dA~ {Ho ' Aj)m + dB~ {Ho ' Bj}m ~ {Ho > fw}m 

(74) 

Hence, in this representation the quantum behavior is displayed by the stargenf unction sector alone. 
However, for this system, we also have (let z — qi,pi Vp2 and \z > be a generic eigenket of z with 
associated eigenvalue z ): 

W (^S(\ Z >< Zo\) = &*>(t,P t ,q,p)( Z ~ Z o) = \(t,4>,A,B)( Z (^ £ B) ~ Z 0)\A=A(t^p)AB=B(tm 

= 5( Z (t, A, B) - z o)\A=A(t,q,p)AB=B(t,q,p) = 5 ( Z ~ Z o) (75) 

where in the third step we used the fact that el^^^ -20 ^ = e jfe O(*A-B)- z o) f or z _ q 1 p 1 y p 2 ^ 
(a simple result that follows from eq.(61)). Hence, the former three fundamental stargenfunctions 
display a classical structure and satisfy -^A^(t,p t ,q,p)( z ~ z o) — 0. We conclude that for this system, 
in this representation, the non-trivial (quantum) behavior is displayed by the stargenfunction z = q<i 
alone. 



6 Conclusions 

We addressed the deformation quantization of the non-relativistic particle in the generally parametrized 
form and provided the complete specification of its physical quantum states in three different phase 
space representations. We proved that in one of these representations the distributional sector 
of the system displays a classical causal structure. This result confirms a similar conclusion that 
has been recently obtained for a generic non-relativistic deparametrized system and reinforces the 
point of view that the De Broglie-Bohm formulation is not the unique possible causal formulation 
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of quantum mechanics. In the De Broglie-Bohm theory HS1 E| the source of quantum behav- 
ior is the quantum potential determining a causal (although not fully classical) dynamics for the 
quasidistribution. Furthermore, the theory also displays a non-trivial quantum correction to the 
momentum stargenfunction. On the other hand, in the "causal covariant formulation" presented 
here the quantum effects have been completely removed from the distributional sector (which now 
displays a fully classical causal structure) and the price to pay was the appearance of some further 
(quantum) corrections on the observable's sector of the theory. 

We finish by recalling the point of view according to which the relation between the classical 
statistical and the quantum mechanical formulations of generally covariant systems should be further 
explored as it may provide key physical insights into some of the conceptual problems displayed by 
the quantum version of these systems. From this perspective the deformation methods seem to be 
especially suited. We proved that (at least for the simplest case of the parametrized non-relativistic 
particle) the formal similarities between classical statistical and phase space quantum mechanics 
carry on intact to the generally covariant context. This close analogy is obviously superior to 
the one displayed by the standard operator formulation, and supports the point of view that the 
deformation methods should be further explored as an alternative, conceptually simpler approach 
to the quantization of generally covariant systems. 
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